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Q\ ■ Abstract 



The unification of gauge and Yukawa coupiings within the minimal supersym- 
metric standard model is studied at the two loop level. We derive an expression for 
the effective scale, Tsusy, which characterizes the supersymmetric particle thresh- 
old corrections to the gauge couplings, and demonstrate that Tsusy is only slightly 
dependent on the squark and slepton masses, and strongly dependent on the Hig- 
gsino masses as well as on the mass ratio of the gauginos of the strong and weak in- 
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5h ' teractions. Thus, if the gaugino masses proceed from a common soft supersymmetry 



C3 



breaking term at the unification scale, and there is no significant gaugino-Higgsino 
mixing, Tsusy will have a strong dependence on the supersymmetric mass param- 
eter \i. Moreover, the value of the top quark Yukawa coupling necessary to achieve 
the unification of bottom and tau Yukawa couplings is also governed by Tsusy , and 
it yields predictions for the top quark mass which are close to the quasi infrared 
fixed point results associated with the triviality bounds on this quantity. From the 
requirement of perturbative consistency of the top quark Yukawa sector of the the- 
ory, we obtain constraints on the allowed splitting of the supersymmetric spectrum, 
which, for certain values of the running bottom quark mass, nib, are stronger than 
those ones coming from the experimental constraints on the strong gauge coupling. 
For example, for m&(M&) < 4.1 GeV, which approximately corresponds to a physi- 
cal mass Mf, < 4.7 GeV, we obtain that, for moderate values of tan/3, perturbative 
unification may be achieved only if a%(Mz) < 0.124. 



*On leave of absence from the Univ. of Warsaw. 
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1 Introduction 



The Standard Model (SM) provides a very good understanding of the strong and 
electroweak interactions and, so far, it has withstood all the experimental onslaugths. 
Yet, there are a host of open questions, which need to be answered. Some of them are 
of technical nature, since they are related to the probation of the model at a higher 
level of refinement upon more precise measurements from the experiments and more 
meticulous theoretical computations to analyze the discrepancies between observations 
and predictions. Apart from this type of questions, there are others which involve more 
fundamental issues, and are those concerning the generation of masses with the associated 
large number of unknown parameters, the properties of matter with the unexplained 
apparent triplication of fermions into families and the peculiar assignment of quantum 
numbers and, most challenging, the mystery of forces. In fact, the SM does not comprise 
all the forces observed in nature. Moreover, even if there were no other forces in nature, let 
aside gravity, besides those associated with SU(3) C x SU(2)i x U(l), it would be natural 
to ask about the physical principle which selects this group structure above others. As a 
consequence, the idea of the existence of an underlying structure unifying all the gauge 
interactions, from which the SM is derived as an effective theory, appears as a sensible 
alternative. In fact, using the Renormalization Group (RG) evolution to determine the 
unknown high energy behavior of the SM parameters, it is possible to envisage that the 
three different gauge couplings may converge to a common value at some unexplored 
high energy scale A 0,0. Indeed, the unification of couplings is the greatest conceptual 
advantage of a Grand Unified Theory (GUT). Other quite important consequence of a 
GUT is that it naturally implies charge quantization. 

The efficiency of the RG approach to explore the short distance physics is, of course, 
strongly dependent on the accuracy of the low energy data to be used as the initial bound- 
ary conditions in the RG evolution. In this respect, considering the available accurately 
determined value of the electromagnetic gauge coupling as well as the precision mea- 
surements for the weak couplings from LEP and from neutrino scattering experiments 
and in spite of the rather large experimental uncertainties in the determination of the 
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strong gauge coupling, it has been already shown in the literature that the simplest GUT 
theories, including SU(5), are ruled out. On the contrary, after the inclusion of super- 
symmetry, the gauge couplings of the SM converge to a common value at a high energy 
scale which may be sufficiently high to prevent a too fast proton decay rate, incompatible 
with experiment 0-0. 

As a matter of fact, recently there has been renewed interest in analyzing models of 



GUT's, and, indeed, some possible high energy scenarios have been studied flqj - |T2 
However, most important is to observe that all these various grand unification scenar- 
ios lead to similar low energy predictions. Thus, for the extent of this work, we shall 
concentrate on the minimal supersymmetric extension of the SM in the light of the mini- 
mal presumption of Grand Unification, without any assumptions about the physics above 
the unification scale. In particular, we shall also analyze the unification of couplings in 
the Yukawa sector, which appears naturally in some grand unified models |l0||-||12|. The 
equality of the bottom quark and r lepton Yukawa couplings is obtained in models with 
SU(5) gauge symmetry, while the extreme case of equality of the r lepton and the top 
and bottom quark Yukawa couplings may occur in SO(10) and E(6) unification schemes. 
Beyond all the details depending on a particular high energy model, which we shall not 
address in the present study, what remains effective is that definite predictions for the 
top quark mass can be made in the framework of gauge and Yukawa coupling unification. 

The issue of the unification of couplings in the Minimal Supsersymmetric Standard 
Model (MSSM) has been addressed to various extensions in previous analyses P|-[|T2" 



However, in the present work, we study some aspects of the Grand Unification framework, 
which yield new interesting results in relation to the supersymmetric particle spectrum 
and the supersymmetric threshold usually considered in the literature. Our paper is 
organized as follows. In section 2 we present the experimental predictions, which we shall 
consider to constrain the values of the low energy parameters of the theory. In section 3, 
we describe the two loop RG method that we use in this study. In particular, we explain 
in detail the treatment of the different threshold corrections which are incorporated in 
our analysis. Most important, we derive an expression for the effective scale, Tsusy, 
which governs the supersymmetric threshold corrections as a function of the different 
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super symmetric particle masses. In section 4 we discuss the general features of gauge 
and Yukawa coupling unification. We present an approximate analytical computation 
for the top quark mass derived from the infrared quasi fixed point associated with the 
triviality bounds on the top quark Yukawa coupling, and compare these results with our 
numerical two loop predictions. In section 5 we analyze the constraints on Tsusy and, 
therefore, on the splitting of the sparticle spectrum, which, depending on the exact values 
of the bottom quark mass and the weak mixing angle sin 2 9yy, may be derived either from 
bounds on the top quark Yukawa coupling coming from the requirement of perturbative 
consistency, or from the experimental bounds on the strong gauge coupling, a^iMz)- 
Section 6 is devoted to the study of the two loop RG results for the top quark mass as a 
function of tan f3 - the ratio of the vacuum expectation values of the two Higgs fields, and 
its dependence on a 3 (M z ) (or equivalently on sin 2 9w) an d the bottom quark mass. In 
particular, we analyze the conditions for which the three Yukawa couplings of the third 
generation acquire a common value at the unification scale, Mqut- We also present the 
predictions for the light Higgs mass as a function of tan/3 and their dependence on the 
squark mass. We reserve section 7 for our conclusions. 

2 Experimental Constraints on the Low Energy Pa- 
rameters 



In order to analyze the unification of gauge couplings in the MSSM, and obtain reliable 
results from the RG analysis, it is necessary to consider the experimental data, to deter- 
mine the low energy parameters as precise as possible. Due to the most recent analysis 
of the LEP data, the electroweak parameters are known at present with high precision. 
The largest uncertainty comes from the unknown top quark mass, which induces not only 
logarithmic but also quadratic corrections on the value of sin 2 #y^. In the following, we 
shall define the gauge couplings in the modified minimal substraction renormalization 



scheme (MS), introduced in Ref.JT3|. In such case, the best fit to the data is obtained 
fori 

sin 2 9 W (M Z ) = 0.2324 ± 0.0003, (1) 
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for a top quark mass value M t = 138 GeV. If the top quark mass is left free, the central 
value does not change, but the errors increase. For other values of the top quark mass, 
the best fit to the data is obtained for values of the weak mixing angle which are related 
to M t by 

sin 2 6 W (M Z ) ~ 0.2324 - 1. 1(T 7 GW~ 2 (m 2 - M 2 ) , (2) 

with M to = 138 GeV. This implies that, for a top quark mass M t ~ 100 GeV, the best 
fit is for sin 2 8 W ~ 0.2335, while, for M t ~ 200 GeV, one obtains sin 2 6 W ~ 0.2305. The 
value of the electromagnetic gauge coupling is only logarithmically dependent on M t and 
for Mt = M to is given by 

' = 127.9. (3) 



a(M z ) 

On the contrary to what happens in the electroweak sector of the theory, the strong 
gauge coupling a^{M z ) is not so precisely known. A general tendency in the experimental 
determination of a^(Mz) is that those values which are obtained from low energy exper- 
iments are slightly lower than those obtained from LEP data. Although a light gluino 



could improve the agreement between the low and the high energy data fll3 |, there are 



theoretical uncertainties in both determinations, and these uncertainties may account for 



the apparent discrepancies in the value of a^iMz) [0- A conservative attitude is to take 



the value for the strong gauge coupling to be in the range 

a 3 (M z ) = 0.12 ±0.01. (4) 

The top quark has not yet been observed, but, as mentioned above, within the SM, its 
mass may be indirectly determined through its radiative correction effects, in particular, 
through the measurement of the p parameter. Conservative upper and lower bounds on 
the physical top quark mass are given by 

91 GeV <M t < 200 GeV. (5) 

The bottom quark, instead, has been already observed experimentally. However, there is 
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still a large uncertainty about its physical mass. Following the particle data book fl6 
shall take values for the physical bottom quark mass in the range, 

M b = 4.7 - 5.2 GeV. (6) 



One should remember, that there is a significant quantitative difference between the 
running bottom (top) quark mass and the corresponding physical mass, defined as the 
location of the pole in its two point function. This difference comes from QCD corrections, 
which, at the level of two loops, are given by, 



where M q and m q are the physical and running quark masses, with q = b,t denoting the 
bottom and top quarks, respectively, and Kb = 12A, K t ~ 11 flT7|. From Eq.(^), we 
obtain that, for a^{Mz) = 0.115 (0.125) and a physical bottom quark mass = 4.7,4.9 
and 5.2 GeV, the corresponding running mass parameter is m^M?,) ~ 4.1(4.0), 4.3(4.2) 
and 4.6 (4.5) GeV, respectively. The bottom quark running mass is typically a 15% lower 
than the physical mass, with a slight dependence on the exact value of the strong gauge 
coupling a 3 (Mz). As we shall discuss below, assuming bottom and tau Yukawa coupling 
unification, the difference between the physical and running bottom quark masses has 
relevant consequences on the top quark Yukawa sector. The reason is, that the bottom 
quark mass value fixes the overall scale of the bottom quark Yukawa coupling. Our 
predictions are, instead, only slightly dependent on a moderate variation of the energy 
scale at which the running mass is evaluated. 

In contrast to the third generation of quarks, the tau lepton mass is very well known, 



In this case, the difference between the running and the physical masses is only due 
to electromagnetic corrections and can be neglected for the aim of the computations 
performed in the present work. 



3 Renormalization Group Analysis and Threshold Cor- 
rections 





M T = 1.78 GeV. 




A first approach to the problem of unification of couplings within the context of the 
MSSM is to assume the absence of threshold corrections at the grand unification scale. In 
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addition, a usual assumption considered in the literature is that all the supersymmetric 
masses are fixed at a common energy scale, M SUSY - In such case, the effective scale Tsusy, 
which characterizes the supersymmetric threshold corrections to the gauge couplings, is 
directly equal to Msusy- Once these assumptions are made, for a given top quark mass, 
the value of the weak mixing angle may be predicted as a function of the strong gauge 
coupling a 3 (M z ) and the overall soft supersymmetric threshold scale Msusy = T S usy- 
Within the above framework, it has been shown that, for values of the weak mixing 
angle and the strong gauge coupling consistent with experimental data, the unification of 
couplings may be achieved for values of the grand unification scale of the order of 10 16 GeV 
and Msusy of the order of the weak scale |§. Based on this fact, it has been first thought 
that a more precise determination of the values of a^{Mz) and sm 2 6w(Mz) would result 
in a prediction for Tsusy and, hence, for the characteristic value of the supersymmetric 
particle mass M S usy- However, further theoretical work 0-0, has shown that, once 
the condition of equality of all supersymmetric particle masses is relaxed, the effective 
supersymmetric threshold scale Tsusy can be very far from any of the sparticle mass 
scales. Thus, the meaning of Tsusy as the parameter which represents the characteristic 
sparticle mass scale is no longer valid. Nevertheless, as we shall explicitly show below, 
the consequences of an arbitrary sparticle spectrum on the predictions for the measurable 
low energy gauge couplings can always be parametrize in terms of a single threshold scale 
Tsusy 0- 

Due to the significance of the effective scale Tsusy for our study, it is relevant to 
analyze its behavior in detail. Langacker and Polonsky [0] have derived a formula which 
determines Tsusy in terms of three low energy effective threshold scales Tj associated 
with each gauge coupling on []■ However, in order to obtain more information from the 
dependence of our predictions on the effective supersymmetric threshold scale, we find 
it most useful to present Tsusy in terms of the physical masses of the supersymmetric 
particles. 

1 Observe that there is a slight difference between our notation and that one used by the authors of 
Ref. §. The M. t scales there are our Tj and they have used the parameter Msusy to denote both the 
meaning of our parameters Msusy and Tsusy- 
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3.1 Threshold Correction Effects 



We assume the unification of the three gauge couplings at Mqut and consider the beta 
functions for their running up to two loops. Thus, after the inclusion of the one loop 
threshold corrections coming from the supersymmetric particles, for each gauge coupling 
ctj, we obtain the relation 



1 1 b fSSM f M GUT \ 1 

a G OLi{M z ) 27r V M z J af tnr 
where ji contains the two loop contributions to the beta function, while 

* V S.in(^), (10) 

are one loop threshold corrections to l/a^M^). In the above, the summation is over all 
sparticle and heavy Higgs doublet states with masses above Mz-, is the contribution 
of each sparticle (or heavy Higgs) to the one loop beta function coefficient of the gauge 
coupling ttj and b^ SSM is the one loop beta function coefficient of the gauge coupling 
oti in the MSSM. Finally, the term Aj includes the threshold corrections related to the 
top quark mass as well as a constant term which depends only on the gauge group, 
Ai = 0, A 2 = -1/67T, A 3 = — 1/4-7T, and results from the necessity of converting the 
gauge couplings, which are evaluated in the MS scheme, to their expressions in the 
Dimensional Reduction (DR) scheme, in order to work consistently in the MSSM fI8 |. It 
is important to remark that, when solving the renormalization group equations at the two 
loop level, it is in general sufficient to consider threshold corrections at the one loop level. 
In the above, we have not included threshold corrections arising at the grand unification 
scale. The quantitative study of these corrections requires the knowledge of the model 
describing the physical interactions at the high energy scale, something that is beyond 
the scope of the present analysis. We shall elaborate on this issue in a future work |19| . 

The unification condition allows to determine one coupling, for example a^Mz), as 
a function of the electroweak gauge couplings at M z and the one loop supersymmetric 
threshold corrections. Following first the parametrization of Langacker and Polonsky 0, 
let us define three low energy effective threshold scales Tj, associated with the one loop 
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super symmetric threshold corrections to the gauge couplings ctj by 

where 6f M is the one loop beta function coefficient of the gauge coupling aij in the SM. 
Observe that, with the above definition, Eq . (|i~l~D , if all sparticles acquire masses above the 
scale Mz the effective scale 7] will be larger than Mz and of the order of the characteristic 
mass of the sparticles contributing to hf ISSM . From Eqs. (||) - (PI); it follows that, 



a 3 (M z ) ' (h-b 2 ) 



+ 72 + A 2 



a 2 (M : 



z) 



(&2 - h) 
(h - b 2 ) 



+ 7i + Ai 



73 -A 3 + A^ — -— (12) 
\a 3 {M z ) J 



where 

A «fcr ( 1 \ _ (fel-&3) 1 _ fa-h, ■ , , , 



a 3 (M z ) ; (6 X - b 2 ) al thr (b, - b 2 ) af thr af hr ' 



is the contribution to \ja 3 (Mz) due to the inclusion of the supersymmetric threshold cor- 
rections at the one loop level. In the above, for simplicity of notation, we have dropped the 
superscription from the supersymmetric beta functions, naming bf ISSM = b{. Replacing 
the values of bi and bf M in Eq.([U|), while using Eqs.(|II]) and (|TT|), we obtain 



A 



Sthr 



a 3 (M z ) J 28tt 



100 In {—) -251nf— ) -56h/ ' 



(14) 



M z ) \M Z ) \M Z , 

Thus, the supersymmetric threshold corrections, Eq. ([T3|) are characterized by the mass 
scales Tj. 

3.2 The Effective SUSY Threshold Scale 

The effective threshold scale Tsusy is defined as that one which would produce the same 
threshold correction given in Eq. flTj]) , in the case in which all the supersymmetric particles 
were degenerate in mass. Hence, as a function of the mass scales Tj, it reads, 

100 In (%L) - 25 In ( £-) - 56 In ( A) = 19 In , (15) 

\M Z ) \M Z ) \M Z ) \ M z J V ; 

in agreement with the expression given by Langacker and Polonsky M. From equation 



(15") we can easily see that, unless 7\ ~ T 2 ~ T 3 , the effective supersymmetric threshold 



scale Tsusy may significantly differ from the actual value of the masses Tj. In particular, 
it may be lower than Mz- 

The low energy input parameters for the RG equations relevant to achieve gauge 
coupling unification are given at the scale Mz- Thus, it is necessary to understand, which 
is the role of Tsusy in the case in which its value occurs to be below Mz- In fact, if, for a 
fixed set of values of the low energy coupling constants ai(Mz), the unification of couplings 
takes place for a given splitting of the sparticle spectrum equivalent to consider a value of 
Tsusy > M z , then, for the same values of ctj(M^), other different splittings of the sparticle 
spectrum will still be compatible with gauge coupling unification, as far as the value of 
the effective threshold scale, as defined in Eq . (|i5|) , is preserved. In particular, the case 
in which all sparticles acquire a mass Msusy = Tsusy will be the possible alternative 
for vanishing splitting of the supersymmetric spectrum. If Tsusy < Mz, instead, for 
a given set of low energy parameters the unification of the gauge couplings at a scale 
Mqut will only be achieved in the case of a sufficiently large sparticle mass splitting, 
which, through Eq.fll^), gives the corresponding low value for T S usy- Moreover, T S usy 
conserves a proper mathematical meaning: Unification of couplings will be achieved if we 
run down the standard model RG equations from Mz to Tsusy and then we run up the 
supersymmetric RG equations from Tsusy to the unification scale. On the contrary, if 
all sparticles had a characteristic mass M S usy lower than M z , they would produce no 
threshold corrections to the values of ai(M z )- Thus, in such case T S usy can not be but 
Mz, and, therefore, it is no longer possible to identify it with Msusy as before. This 
means that, for values of the low energy couplings, for which we obtain unification for a 
splitted supersymmetric spectrum which yields a value of Tsusy < Mz, we would never 
obtain unification in the case of a degenerate supersymmetric spectrum. 

A remark is in order. The scale Tsusy is only related to the condition of unification 
and not to the specific value of M GU t and o.c{Mqut) ■ In fact, for the same value of T S usy, 
different sparticle spectra will produce the unification of the gauge couplings at different 
grand unification scales and with different values of the coupling constant aa(MauT)- 
This has implications for other constraints, coming, for example, from the nonobservation 
of proton decay [§] - |§ . 
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3.3 Tsusy and the Physical Sparticle Masses 



28/19 

/ / / / . ; \ 

TsUSY = TTlg 



m w 



m-g 



In order to study the dependence of Tsusy on the different sparticle mass scales of the 
theory, we define m^, trig, m z ~, m^, and ran as the characteristic masses of the 
squarks, gluinos, sleptons, electroweak gauginos, Higgsinos and the heavy Higgs doublet, 
respectively. Assuming different values for all these mass scales, we derive an expression 
for the effective super symmetric threshold Tsusy which is given by 

\ 3/19 / \ 3/19 / \ 4/19" 

s H N . (i6) 

™§J \ m Hj \ m Hj 

The above relation holds whenever all particles considered above have a mass m, n > Mz- 
If, instead, any of the sparticles or the heavy Higgs boson has a mass < Mz, it should 
be replaced by Mz for the purpose of computing the supersymmetric threshold corrections 
to \/az{Mz)- In the following, unless otherwise specified, we shall assume that all sparti- 
cles and the heavy Higgs doublet acquire masses above Mz- From Eq . (|16"D , it follows that, 
for fixed mass values of the uncolored sparticles, that is sleptons, Higgsinos and the weak 
gauginos, together with the heavy Higgs doublet, the value of T S usy decreases for larger 
mass values of the colored sparticles - squarks and gluinos. Moreover, Tsusy depends 
strongly on the first two factors in Eq.flTB]), while it is only slightly dependent on the ex- 
pression inside the squared brackets. This is most surprising, since it implies that Tsusy 
has only a slight dependence on the squark, slepton and heavy Higgs masses and a very 
strong dependence on the overall Higgsino mass, as well as on the ratio of masses of the 
gauginos associated with the electroweak and strong interactions. The mild dependence 
of the supersymmetric threshold corrections on the squark and slepton mass scales is in 
agreement with a similar observation made in the context of the minimal supersymmetric 
SU(5) model in Ref.g. 

In the above, we have implicitly identified the sparticle mass eigenstates with SU (3) c x 
SU(2)l x U(l)y eigenstates. Since the most significant dependence of Tsusy on the 
supersymmetric particle masses comes from the gaugino and Higgsino sectors, we expect 
no significant dependence of Tsusy on the squark and slepton mixing. The neutralino 
and Higgsino mixing may, instead, have quantitative effects, which have not been included 
above. We shall deepen into this point elsewhere ||19| . 
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In some models, in which the source of supersymmetry breaking is given by a common 
gaugino mass, M1/2, at high energy scales, it follows that, within a good approximation, 
all colored sparticles have a common mass m co \ = 0(rrig) while all uncolored sparticles 
acquire a common mass m W k = 0(Mi/ 2 ). The mass of the Higgsinos is determined by the 
super symmetric mass parameter /i appearing in the superpotential and we are implicitly 
assuming that \fi\ = 0(Mi/ 2 )- Thus, in the simplified case in which only two different 
supersymmetric mass scales are presumed, the effective supersymmetric threshold scale 
is given as a function of both mass parameters m co i and m W k as follows, 



As we remarked above for the most general case, if any one of these two masses is below 
M z , it should be identified with M z while applying Eq. ([17|) . Assuming m co i > m wk , 
Tsusy is predicted to be smaller than m wk . Moreover, since Tsusy decreases for larger 
values of m co i we have that for a fixed value of m w k, the larger the splitting of masses the 
smaller becomes Tsusy- For example, if all the colored sparticles are assumed to have 
masses of the order of m co i ~ 1 TeV, while the uncolored ones have masses of the order 
of m w k = 200 GeV, the resulting effective supersymmetric threshold scale is as low as 
Tsusy — 15 GeV. Of course, the fact of considering all colored sparticles to be degenerate 
in mass and the same for all the uncolored ones is an approximation, which becomes 
natural only in some particular supersymmetry breaking schemes. However, due to the 
slight dependence of Tsusy on the squark, slepton and heavy Higgs masses, the effective 
supersymmetric threshold scale given in Eq.(jI5|) is approximately described by 



Therefore, the approximation made above, Eq. fll7D , turns out to be a rather good one if 
one identifies m co ; with the gluino mass, while m w k is a geometrical average of the overall 
value irtfj of the Higgsinos and that one of the electroweak gauginos. This simplification 
yields a clearer computational framework for all quantities with a dominant dependence 
on Tsusy- 

Depending on the SUSY breaking scheme, it may occur that there is no mass splitting 
between squarks and sleptons. If, for example, the dominant source of soft supersymme- 







12 



try breakdown is a common mass tuq 3> M§ for the spin-0 sparticles, while the common 
gaugino mass Mi/ 2 is of the order of M z , then squarks, sleptons and the heavy Higgs 
particles will become heavy and almost degenerate in mass, while gauginos will be rel- 



atively light p0 |. However, in such framework as well as in the case in which squarks 
and sleptons are far from being degenerate in mass, the behavior of the supersymmetric 
threshold scale as a function of the sparticle masses implies that Eq. ( p~8|) gives, to a good 
approximation, the value of T S usy- 

In many grand unified scenarios, the supersymmetric spectrum is determined by the 
soft super symmetry breaking parameters arising at Mqut, that is a common scalar mass 
m , a common gaugino mass M1/2, and the scalar trilinear couplings A Q and B involving 
contributions proportional to the trilinear and bilinear terms in the superpotential, as well 
as the supersymmetric mass parameter fiQ, appearing in the superpotential. Moreover the 
assumption of small gaugino - Higgsino mixing implies that the overall Higgsino mass is 
approximately given by the absolute value of the renormalized supersymmetric mass fi, 
while, the assumption of an overall common gaugino mass at Mqut implies that the ratio 
of the gluino to the wino mass is approximately given by the ratio of the strong to the 
weak gauge coupling. Therefore, if all sparticle masses are above Mz, we obtain that the 
approximate expression of Tsusy, Eq.(p!8[) simplifies to 

Thus, the effective supersymmetric threshold scale is proportional to the absolute value 
of the supersymmetric mass fi {Tsusy — 1/^1/7), and it is by no means a measure of the 
scale of soft supersymmetry breaking. Taking values for | /x | < 1 TeV, we obtain that the 
effective threshold scale is of the order of, or smaller than Mz- Observe that, since within 
this framework the threshold corrections to 0:3 {Mz) depend on a^{Mz) itself, the value 
of the strong gauge coupling must be found by solving Eq . (|12D in a self consistent way. 
Within the minimal SU(5) supersymmetric model, this was performed in Ref.jQ. 

Once more, if any of the sparticles acquires a mass which is small in comparison with 
Mz, this mass must be replaced by Mz when applying Eq. (|T8D . For example, if the 
supersymmetry breaking parameter M1/2 were very small in comparison with Mz, all 
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gaugino masses would be smaller than or of the order of Mz- In this case, the effective 
supersymmetric threshold scale would be approximately given by the Higgsino mass. This 
situation is, however, not favored by experiments, apart from the possible existence of a 



light gluino window Ul~4" |. Another possibility is that the winos acquire masses below Mz, 
while the gluinos are heavier than the Z° vector neutral boson. In such case, one should 
replace the wino mass by Mz directly in Eq.(|T8|). 

3.4 Numerical Computations 

Once the low energy values of the fine structure constant, a(M z ), and sin 2 W (M z ) are 
fixed, the condition of gauge coupling unification allows to determine the value of a^iMz), 
up to threshold corrections. For the aim of performing the numerical computations, rather 
than including the supersymmetric threshold corrections through the procedure explained 
above, we shall approach the issue in the following equivalent way: At energies above the 
heaviest sparticle mass, we use the complete two loop renormalization group equations 
of the MSSM, including the effect of the SU(3) C x SU(2) L x U{1) Y gauge couplings and 
the third generation of quark and lepton Yukawa couplings. Below the scale of energy 
given by each sparticle mass, we account for the threshold corrections by modifying the 
one loop renormalization group equations while keeping the two loop contributions as in 
the complete MSSM (this is a correct procedure, up to higher order corrections). We 
shall assume, as in Eq . (|T7|) , that all colored sparticles acquire a common mass m co i, while 
the uncolored ones, together with the heavy Higgs doublet acquire a mass m wk . Thus, 
in practice, we consider a two step decoupling for the supersymmetric particles and the 
heaviest Higgs doublet. At energy scales m w k < \i < m co i, the model has the same particle 
content as in the standard model with two Higgs doublets and with the addition of all 
the uncolored sparticles. At energies below m„i, we recover the standard model particle 
spectrum. Thus, we account for the uncolored sparticle and heavy Higgs doublet threshold 
corrections while decoupling them at m W k and using the full two loop renormalization 
group equations of the SM for the energy range Mz < < m w k. Furthermore, since the 
top quark mass fulfills the relation M t > Mz, the top quark threshold corrections must 
be also taken into account. We include the top quark threshold corrections to the gauge 
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couplings analytically, by the procedure explained in Ref.0. 

The inclusion of the Yukawa sector in our study requires extra boundary conditions for 
our numerical analysis, besides the low energy electroweak parameters and the condition 
of gauge coupling unification already implemented in the gauge sector. Indeed, the values 
of the bottom quark and tau masses are necessary input data for the low energy bound- 
ary conditions while the unification of their Yukawa couplings determines the boundary 
conditions at the unification scale. Then, we obtain predictions for the top quark mass as 
a function of tan f3. In the extreme case in which we require the unification of the Yukawa 
couplings of the top and bottom quarks and the tau lepton at the unification scale - 
h t (M GUT ) = h b (M GUT ) = h T (M GUT ) - a prediction for tan/? is derived as well. Obviously, 
while performing the RG evolution in the Yukawa sector, in order to be consistent within 
our framework, we have to consider the same two step decoupling procedure at m co i and 
m wk and to include the threshold corrections associated with the decoupling of the top 
quark, too. At energy scales /x < M z , we have an effective SU(3) C x U(l) em symmetric 
theory. Hence, we consider the running of the bottom quark and r lepton masses within 
the SU(3) C x U(l) em gauge group. 

3.5 The Higgs Sector 

We also consider the Higgs sector of the theory |21| - [29]|. At energy scales above m co i, 



the Higgs potential in the MSSM with general soft-SUSY breaking terms is given by |2 



29 



V eS = m\E[H x + m\E\E 2 - mj(HfiT 2 H 2 + h.c.) 

+ gZ)(B$H 3 - H\Hi) 2 + isflfljffxl 2 , (20) 

where g 2 and g\ are the gauge couplings of SU(2)l and U(l), respectively. The mass 
parameters mj, with i — 1,2, 3, as well as the gauge couplings are running parameters 
depending on the renormalization scale. The parameter M 2 , defined as m 2 + ml = 2M 2 , 
must be of the order of m w k, to be compatible with the assumption that only one light 
Higgs doublet appears in the low energy spectrum of the theory. At energy scales /i < m coh 
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the potential is given by the general expression 

Ves = mlHlH 1 + mlHlH2-ml(H^iT 2 H2 + h.c.) 

+ ^ (HlH,) 2 + ^{H\H 2 f + X^HlH^HlHz) + \,\H\it 2 HI\ 2 . (21) 

Therefore, the running quartic couplings Xj, with j = 1 — 4, must satisfy the following 
boundary conditions at m co f. 

a 9i + 92 \_92~9i \ 9 2 / 00 \ 

Ai — A 2 — , A 3 - , A 4 - — — , {11) 

where we have used the relation 

{H\H x )(HlH 2 ) = \H\Hi\ 2 + \H\iT 2 Hl\ 2 . (23) 

At energy scales /i below m coh the values of the quartic couplings may be obtained by 
solving the corresponding renormalization group equations |29H,|l30|, considering the step 



decoupling described above. As we explained in Ref . |26| , the appearence of other quartic 



couplings than those ones given in Eq.(|21|) is protected by either discrete symmetries 
(H 2 — > — H 2 ; Hi — > —Hi) or a global PQ symmetry. In particular, a quartic term 
X^KH^^Hi) 2 + h.c] breaks the PQ symmetry. The PQ symmetry is only broken by 
mass or scalar trilinear terms, for example the term m 2 appearing in the scalar potential, 
Eq. (|2T|). Hence, A5 will not receive any leading logarithmic contribution and will not 
appear in the renormalization group analysis. 

Quite generally, after the electroweak symmetry breaking, three neutral and two 
charged scalar states appear in the spectrum. As we said above, in our analysis we 
are considering the case in which only one light Higgs doublet $, which is a combination 
of the original Higgs doublets Hi and H 2 , remains below m^. It is straightforward to 
prove that the mass matrices of the charged and CP-odd states are always diagonalized 
with an Hi — H 2 mixing angle which is given by —f3. In addition, whenever the mass 
parameter M 3> Mz, the mixing angle for the neutral CP-even states is approximately 
given by /3 and, hence, the lightest CP-even state, together with the Goldstone modes 
form a Higgs doublet whose expression is given by 

$ = Hi cos p + ir 2 H* sin (3. (24) 
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The mass of the lightest CP even Higgs particle is given by 

m\ = 2A(m> 2 (25) 

where v 2 = v \ + t>f, with t> 12 the vacuum expectation values of the Higgs fields H 12 , 
respectively, and A(/x) is the effective quartic coupling of the light Higgs doublet. The low 
energy value of the quartic coupling A may be obtained by using the RG equations with 
the boundary condition 

A = Ai cos 4 /3 + A 2 sin 4 (3 + 2(A 3 + A 4 ) cos 2 (3 sin 2 (3 (26) 

at the scale \x = m w k- In Eq.(^) we have defined the light Higgs running mass at m t , 
which, in a very good approximation coincides with the on shell mass m/^m/J. This 
is due to the fact that the only relevant radiative corrections to the Higgs particle are 
induced by the top quark, and in the MSSM, for values of the top quark mass compatible 
with experimental bounds, m t > Mz, the relation m t > rrih is fulfilled. Therefore, the 
evaluation of the light Higgs running mass at mt allows a correct evaluation of the top 
quark threshold corrections, while, below m t the Higgs mass running is negligible. 

The full two loop RG equations we use in the MSSM and the SM regimes, are given 
in Ref.plf and Refs.[|17|, P2| , respectively. In addition, the RG equations for the bottom 



quark and tau lepton considered at scales /x < M Z} are given, for example, in Ref . flTf . The 



proper one loop renormalization group equations, describing the evolution of the gauge, 
Yukawa and quartic couplings, in the case in which the uncolored sparticles are the only 
supersymmetric contributions to the particle spectrum, as it occurs at the intermediate 
scales m co i < \i < m w k, are given in Refs. |2l| , [0. We solve the RG equations numerically 
using the above described multiple step decoupling procedure. 



4 Unification of Couplings 

As we mentioned in the introduction, there have been several recent studies dealing with 
the problem of gauge and Yukawa coupling unification within the MSSM. One of the most 
interesting results is that the values of the gauge couplings obtained from the most recent 
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experimental data is consistent with the unification of gauge couplings within the MSSM, 
if the supersymmetric threshold scale is of the order of the weak scale. Although, quite 
generally, this statement is correct, there are several important points which should be 
taken into account. In what follows, we shall assume that the unification of the weak and 
strong gauge couplings takes place, and we shall leave the inclusion of threshold correc- 
tions arising at the GUT scale for a future analysis |19[. In such context, the following 
properties are fulfilled: 

I) For a given value of T S usy, the predicted value of a 3 (M z ) depends strongly on the exact 
value of sin 2 9w(M z ), and increases (decreases) when sin 2 9w(Mz) decreases (increases). 

II) As can be seen from Eqs.(|l2l)-(|i"5|), for a given value of sin 2 9w(Mz), the predicted 
value of as(Mz) is only logarithmically dependent on Tsusy and increases (decreases) 
when T S usy decreases (increases). 

III) At the two loop level, the contributions of the top quark Yukawa coupling to a 3 (M z ) 
become relevant. For example, for the central values sin 2 8 ] w(Mz) = 0.2324, M t = 138 
GeV, Tsusy — Mz and tan/? such that the top quark Yukawa coupling h t (m t ) ~ 
ht(McuT) ~ 1, the predicted value of the strong gauge coupling reads a^{Mz) — 0.124 
0. However, since at the two loop level the strong gauge coupling receives negative con- 
tributions from the top quark Yukawa coupling, a larger value of the top quark Yukawa 
coupling at Mqut will imply a lower value for a 3 (M^). In fact, a%(Mz) may differ in up 
to a 2% while varying the top quark Yukawa coupling from htiMcvT) — 1 to very large 
values, h 2 t {M GUT ) / '4vr » 1. 

To analyze the consequences of properties I)-III), one must remember that, as we 
discussed in section 3, the effective scale Tsusy is directly related to the overall scale 
of sparticles masses only in the unlikely case of an absolute degenerate supersymmetric 
spectrum. If a significant mass splitting is allowed, Tsusy can significantly differ from the 
characteristic sparticle mass scale. In particular, even when all sparticle masses are larger 
than Mz, the effective Tsusy may be well below the Z° mass scale. Indeed, for a given 
value of sin 2 9w(M z ), property II) implies that it is possible to obtain a lower bound on 
Tsusy by requiring a 3 (M^) to be below its present experimental upper bound 

a 3 (M z ) < 0.13. (27) 
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In principle, one could obtain an upper bound on Tsusy from the experimental lower 
bound on the strong gauge coupling (a 3 (M z ) > 0.11), as well. However, in the present 
analysis, in which no threshold corrections at the unification scale are included, the value 
of Tsusy needed to obtain values of a^{Mz) below its experimental value is, in general, 
of the order of or greater than 1 TeV and, hence, it does not yield significant constraints 
on the sparticle spectrum, once the absence of fine tunning is required. For example, 
for sin 2 6 W {M Z ) = 0.2335 and T SUSY = 1 TeV, we obtain a 3 {M z ) ~ 0.111. Thus, the 
constraint Tsusy — 0(1 TeV) follows in this case. For smaller values of sin 2 QyjiMz) the 
constraints on Tsusy are even weaker. 

For given low energy values of the gauge couplings, the requirement of bottom and tau 
Yukawa coupling unification determines the values of the top quark mass as a function of 
tan/?, depending on the input value of the bottom quark mass |l0|- ||12|| . It has been ob- 



served |]lC)l|,[p6] that, for running bottom quark masses m b (M b ) lower or equal to 4.6 GeV, 
the top quark mass obtained in these analyses is remarkably close to its infrared quasi 



fixed point value |33j associated with the "triviality" bound on the top quark Yukawa 
coupling. As a matter of fact, once the QCD corrections are taken into account, from 
Eq. (0) , it is clear that for values of the bottom quark mass consistent with the experimen- 
tal observations, M b < 5.2 GeV, the running bottom mass is bounded to be rrib(Mb) < 4.6 
GeV. Therefore, the top quark Yukawa coupling obtained in the context of bottom and 
tau Yukawa couplings unification is, indeed, remarkably close to its infrared quasi fixed 
point value. Moreover, a value of the top quark Yukawa coupling close to its infrared 
quasi fixed point necessarily implies that it is getting large at scales of the order of the 
grand unification scale |26| . 



4.1 Infrared Quasi Fixed Point Predictions 

The results associated with the infrared (IR) quasi fixed point can be readily obtained 
analytically. As a matter of fact, recalling the expressions for the beta functions of the top 
and bottom quark Yukawa couplings, and neglecting the small contributions of the weak 
couplings and the tau lepton in this approximation, we observe that the existence of the 
IR quasi fixed point pB]|,||33|| implies straightforward relations among a 3 and Y t = h 2 /(4ir) 
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and Y b = hi/ (An), at an energy scale of order m t . Namely, from the one loop beta function 
of the top quark Yukawa coupling we obtain, 



(16/3)a 3 (mt) - 6Yt(m t ) + Y b {m t ). (28) 

Moreover, the bottom quark Yukawa coupling at the top quark mass scale is approximately 
given by 

h b (m t ) ~ \fl + tan/? 2 m b {m t )/v. (29) 

Since m b {m t ) ~ 3 GeV, for moderate values of tan/3, it follows that Y b {m t ) <C a^{mt) 
and, hence, the bottom quark Yukawa coupling contribution to Eq.(^) can be safely 
neglected. This implies that, for small and moderate values of tan/3 < 30, using Eq . (|28"D 
we can determine the value of the top quark Yukawa coupling as a function of the strong 
gauge coupling constant. In addition, using the relation 

ftan/3 

m t {m t ) ~ ht{m t )—j= , (30) 

1 + tan 2 /? 



for a given value of a^{m t ) we obtain the value of the top quark mass only as a function 
of tan/3. For larger values of tan/3 > 30, instead, the bottom quark Yukawa coupling 



contribution to Eq.(28) can no longer be neglected. In fact, for a fixed value of a^(Mz) 
the value of the top quark Yukawa coupling decreases with increasing tan/? till, for a 
sufficiently large value of tan /3, it becomes equal to the bottom quark one. Then, the IR 
quasi fixed point relation reads, 

(16/3)a 3 K) ^ 7Y t (m t ) = 7Y b (m t ). (31) 

For a given value of the strong gauge coupling and the bottom quark mass, we can 
determine the value of tan f3 and hence of the top quark mass, for which the above relation 
is fulfilled. Finally, for even larger values of tan/3, we are in the regime in which Y b ^> Y t , 
and the IR quasi fixed point relation 

(16/3)a 3 (M z )~6Y b (M z ), (32) 

holds. Since, for a cutoff scale A ~ (9(10 16 ) GeV, the quasi infrared fixed point is associ- 
ated with the renormalization group trayectories corresponding to the triviality bounds on 
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the Yukawa couplings, from Eq. (|32]) we obtain a maximun allowed value for tan f3 , above 
which, a perturbative treatment of the bottom quark Yukawa sector becomes inconsistent. 

In Tables l.a - l.c, we present the analytic IR quasi fixed point predictions together 
with those obtained by requiring the unification of the bottom and tau Yukawa couplings 
within our two loop RG analysis, considering sin 2 9w(Mz) = 0.2324, a 3 (M^) ~ 0.122 and 
three different values of the running bottom quark mass : m&(M&) ~ 4.6, 4.3 and 4.1 GeV. 
For these values of the input parameters we compare the results of both methods for the 
values of the top quark mass for three different values of tan f3 as well as for the value of 
tan f3 and the top quark mass at which the three Yukawa couplings of the third generation 
unify. From the above comparison it follows that, for moderate values of tan/3, the top 
quark mass evaluated through the numerical RG approach is only two to four percent 
lower than its value derived from the analytic approximation performed in the context of 
the IR quasi fixed point. Indeed, the values of m t get closer to the infrared quasi fixed 
point predictions (IRP) for lower values of the running bottom quark mass. As we said 
before, the closeness of the top quark mass predictions is related to the relatively large 
values of the top quark Yukawa coupling h t at the grand unification scale. For very large 
values of tan/5 > 40 and larger values of rrib(Mb), the departure of the RG results from 
the IRP is larger, reflecting the lower values of the top quark Yukawa coupling at the 
grand unification scale. For example, for the three values of mft(M&) = 4.6,4.3 and 4.1 
GeV, the unification of the three Yukawa couplings is achieved for Y^Mqut) — 0.04, 0.13 
and 0.43 respectively. Observe that, in the first case, the Yukawa couplings acquire values 
which are of the order of the gauge coupling value q>g(Mgut) — 0.04. 

4.2 On the Perturbative Consistency of the Top Quark Yukawa 
Sector 

For a consistent perturbative treatment of the theory, the top quark Yukawa coupling 
must fulfill the condition /i 2 (/i)/47r < 1 in the whole energy range in which the MSSM is 
considered as a valid effective theory. However, due to the existence of the IR quasi fixed 
point and the behavior of the top quark Yukawa coupling for intermediate energy scales, 
the requirement of perturbative consistency is naturally fulfilled in the low to intermediate 
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energy regime. Thus, what has to be assured is the correct behavior of the top quark 
Yukawa coupling in the high energy regime. Since we are assuming that the MSSM is valid 
up to scales of the order of Mqut, an upper bound, the so called triviality bound on the 
above quantity may be defined as its low energy value consistent with the renormalization 
group trajectory for which it becomes strong at scales of order Mqut, 



h 2 t (M ( 



GUT) 



< 1. (33) 



Air 

This is equivalent to require that the Landau singularity associated with this coupling 
occurs at an energy scale fJ, S i n g > M GUT . As we shall analyze in detail in section 6, the 
triviality bound on the top quark Yukawa coupling h t determines an upper bound on m t 
as a function of tan (3, which, for values of Mqut ~ 10 16 GeV, is approximately given by 
the infrared quasi fixed predictions for this quantity. 

As a matter of fact, since the value of the top quark Yukawa coupling obtained within 
the framework of bottom and tau Yukawa coupling unification may become quite large, 
the requirement of perturbative consistency of the Yukawa sector of the theory, Eq. (|33|) , 
has important implications in defining bounds on the effective supersymmetric threshold 
scale Tsusy- To understand this, we should emphasize that the properties mentioned 
at the beginning of section 4, implying that the value of a 3 (M^) increases (decreases) 
for smaller (larger) values of sin 2 9yy as well as for smaller (larger) values of Tsusy, are 
independent of the top quark Yukawa sector and, hence, are also fulfilled after the inclusion 
of the Yukawa coupling unification condition. In table 2, we present a particular set of 
predictions to support this general result. 

A second point to be stressed is that the interrelation between the top quark Yukawa 
coupling and is not only a property of the infrared, as shown in section 4.1, but, 
there is a strong correlation between Y t (M GUT ) and a 3 (M z ) as well. Our analysis shows 
that, for values of the squark masses, which do not involve a large fine tuning, trig < 3 
TeV, the actual value of Y t (McuT) strongly depends on the value of a^Mz). Indeed, 
once the values of the running bottom quark mass and sin 2 6w are fixed in each running, 
the value of the top quark Yukawa coupling at M GUT depends only on T SUSY and on 
a 3 (M z ). More specifically, Y t (M GUT ) increases (decreases) with a 3 (M z ), and, therefore, 
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it increases (decreases) for smaller (larger) values of Tsusy together with a 3 (Mz). We 
illustrate this behavior in table 3. Moreover, for a given Tsusy, the value of Y^Mqut) 
is almost independent of the spectific values of the supersymmetric particle masses. We 
illustrate this behavior in Table 4, where we fix the value of Tsusy = 1 TeV and show the 
slight variation of Y^Mgut), ois(Mz) and m t (tan/3 = 4) for m co i varying from 1 to 10 
TeV. As readily seen from Table 4, the radiative corrections to the Higgs particle mass, 
instead, depend not only on the top quark Yukawa sector, and hence on T S usy, but they 
have a logarithmic dependence on the squark mass as well. Most interesting is the fact 
that, due to the behavior of a^{Mz) and Y^Mqut) with the effective supersymmetric 
threshold scale, a lower bound on it, Tg^gy, may be derived from the upper bound either 
on a 3 (Mx), Eq . (]27|) , or on Y^Mqut), Eq.(|33|). A detailed analysis of these constraints is 
presented in section 5. 

5 Constraints on the Splitting of the Supersymmet- 
ric Particle Spectrum 

In this section we shall perform a quantitative analysis of the constraints on the super- 
symmetric threshold scale, which may be obtained from the experimental constraints on 
as(Mz), Eq. (p7|) , or from the requirement of perturbative consistency of the top quark 
Yukawa sector of the theory, Eq.Q3"5|). These constraints have a strong dependence on 
the exact value of the low energy parameters of the theory, in particular, of the bot- 
tom quark mass Mb and sin 2 9w(Mz). Concerning the bottom quark sector, we study 
the cases rrib(Mb) = 4.1,4.3,4.6 GeV, which, considering the proper two loop corrections 
for the central value of a^iMz) — 0.12, approximately correspond to a physical mass 
Mb — 4.7,4.9,5.2 GeV, respectively. The RG equations depend dominantly on the run- 
ning mass, having a very slight dependence on the physical mass through the point at 
which the running mass is evaluated. Hence, for larger (smaller) values of a^[Mz)i our 
results will correspond to slightly larger (smaller) values of the physical mass Mb-, which 
are easily obtainable from Eq. (]?]). 

Although highly more accurate than the bottom quark mass prediction, the experi- 
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mental prediction for sin 2 9w(Mz), Eq.(H), depends quadratically on the top quark mass. 
A way to deal with this difficulty is to take the value of sin 2 Ow(Mz), which provides the 
best fit to the experimental data with a free top quark mass value. Consequently, we have 
considered that the weak mixing angle takes values in the range 

sin 2 9 W (M Z ) = 0.2315 - 0.2335, (34) 

which approximately corresponds to the range of values obtained by such fit at the two 
sigma level J7[. Even though we shall not use the relation given in Eq.(Q) to correlate the 
value of the top quark mass and the weak mixing angle, we should mention that the best 
experimental fit to the top quark mass corresponding to the above values of the weak 
mixing angle are M t ~ 170 — 100 GeV, respectively. 

To perform our analysis we proceed as follows: For fixed values of a(M z ), sin 2 Q w (Mz), 
■mb(Mt)) and m T (M T ) and a given effective supersymmetric threshold scale Tsusy, we 
require the unification of the gauge couplings and the bottom quark and tau Yukawa 
couplings to determine a^iMz) arid the top quark Yukawa coupling h t (m t ), respectively. 
We use the latter, h t (m t ), to predict the on shell value of the top quark running mass 
m t (m t ) as a function of tan/3. It is worth remarking that, for moderate values of tan/3 < 
30, the value of the top quark Yukawa coupling obtained from the unification condition 
has only a slight dependence on this quantity and, hence, the variation of the top quark 
mass with tan/? arises almost exclusively from the variation of the vacuum expectation 
value of the scalar Higgs H 2 , which couples to the top quark. In this section, we shall 
implicitly assume that the condition tan (3 < 30 is fulfilled. 

Using the renormalization group equations, we vary the scale Tsusy until one of the 
two conditions, Eq. (]27|) or Eq.([33|) is fulfilled. This determines a lower bound on T S usy 
below which, gauge and Yukawa coupling unification lead to either a loss of perturbative 
consistency of the top quark Yukawa sector, or to values of az{M z ) which are too large to 
be consistent with the present experimental constraints on this quantity. In Tables 5. a - 
5.c we present the constraints on the value of the effective supersymmetric threshold scale, 
Tsusy i as a function of sin 2 9w(Mz) for different values of the running bottom quark mass 
m&(M&). We also give the values of a 3 (M z ) and Y t (M GU T) which are obtained for the 
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different set of parameters, to show where the constraints are coming from in each case. 
In particular, we observe that for the mean values under consideration, sin 2 9w = 0.2324 
and mb(Mb) = 4.3 GeV, the upper experimental limit on az{Mz) is obtained together 
with the upper bound on the top quark Yukawa coupling coming from the perturbative 
consistency of the theory, Y^Mgut) < 1- 

Moreover, from the numerical results, fixing the value of rrib(Mb) ~ 4.3 GeV and 
varying sin 2 Ow(Mz), it follows that a larger (smaller) value of sin 2 9 'w(Mz) implies that 
the lower bound on Tsusy comes from a 3 (M^) (Y^Mqut)). From Tables 5.a-5.c, it may 
be noticed that, in general, for larger values of m b (M b ) the bounds on T S usy tend to 
come from a 3 (M^), while for smaller values of m b {M b ) the bound tends to come from 
Yt(McuT)- Concerning the value of the lower bound on Tsusy, it happens that, for any 
fixed value of the bottom quark mass, a larger (smaller) value of sin 2 Qw(Mz), yields a 
smaller (larger) value of T^gy. Indeed, such variation of T^^y with sin 2 9 W (M Z ) may 
be easily understood. We know that a^(Mz) and Y^Mqut) increase for smaller values of 
Tsusy- Moreover, a larger value of sin 2 9 w(Mz) implies a smaller value of a^{Mz) and 
Y^Mgut)- Therefore, if the value of sin 2 6w(Mz) is increased then Tsusy has to decrease 
to render 0:3 {M z ) and Y^Mgut) larger so that one of these quantities reaches its upper 
bound. 

It is also quite interesting to observe that, for smaller values of the bottom quark mass, 
mb(Mfe) ~ 4.1, which correspond to physical masses M& ~ 4.7 — 4.8 GeV, the condition 
of unification of couplings implies a theoretical upper bound on a 3 (M z ), which turns out 
to be stronger than the experimental one. In fact, for all the experimentally acceptable 
values of sin 2 Q\y(Mz), if rrib(Mb) ~ 4.1 GeV perturbative unification of gauge and Yukawa 
couplings may only be achieved if a^(Mz) < 0.124. 

As we noticed in section 4, from Eq. (|i~7D we obtain that, for fixed values of the un- 
colored sparticle masses, Tsusy decreases for heavier colored particles. Hence, for a fixed 
value of the uncolored sparticle mass, a lower bound on Tsusy allows to put an upper 
bound on the values of the mass of the colored sparticles compatible with gauge and bot- 
tom - tau Yukawa coupling unification. In figures l.a-c we plot the lower bounds on Tsusy 
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in the m co i - m w k parameter space, for fixed mb(Mb) and different values of sin 2 Qw(Mz)- 
Observe that, the region to the left of any of the plotted curves is associated with a value 
of Tsusy lower than the allowed minimum value Tg^g Y f° r t ne given set of parameters, 
and it is therefore ruled out. This implies that, for a given value of the running bottom 
quark mass m&(M&) and sin 2 9w{Mz), fixing a value for the uncolored sparticle mass the 
allowed values for the colored sparticle mass are those lying to the right of the correspond- 
ing curve and the value on the curve gives the upper bound on m co i. The sharp change 
to a constant value either for m co i of for corresponds to the region in which one of 
these mass scales is below Mz and is hence replaced by Mz in the expression of Tsusy, 
Eq.( f[7|) . Notice that, the lower the value of Tg^ SY , the less constrained the m co i — m w ^ 
parameter space is. Thus, from our explanation above this implies that the stronger con- 
straints on the values of the colored and uncolored mass spectrum compatible with gauge 
as well as bottom - tau Yukawa coupling unification are obtained for the smaller values 
of sin 2 9 W (M Z ) and the bottom quark mass considered. 

In section 4 we comment on the fact that, at the two loop level the gauge couplings 
receive contributions from the Yukawa sector. In particular, the larger the value of the 
top quark Yukawa coupling the larger is its negative contribution to a^iMz)- In fact, for 
values of Y t (M GUT ) = h^(M GUT )/4:Tc ~ 0.1, 0.2, 0.4, 1.0, the two loop top quark Yukawa 
correction to a^(Mz) is approximately given by A Yuk (a^(Mz)) — —0.001, -0.0015, - 
0.002, -0.003 respectively. Hence, for larger values of Y t (M GUT ), a slightly lower value of 
Tsusy is needed in order to achieve unification for a given fixed value of a 3 (M z ). The 
constraints on Tsusy given in table 5 are based on the condition of gauge and Yukawa 
coupling unification, which, quite generally, for values of m&(M&) < 4.6 GeV leads to rather 
large values of the top quark Yukawa coupling. If we relaxed the condition of Yukawa 
coupling unification then it would be possible to consider smaller values for Y t (M G uT) and, 
hence, the corresponding values of (Mz) would be slightly larger. Therefore, if smaller 
values of Y t (M GU T) were taken, the lower bound on T S usy coming from the experimental 
constraints on a 3 (M^) would become stronger due to the increase of the strong gauge 
coupling, in an amount which can be easily worked out by using A Yuk (as(Mz)) and Eqs. 



12|) - (P~5|) . For instance, the value of the strong gauge coupling obtained for sin 6 W (M Z ) 
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0.2324 and T SU sy = 90 GeV at Y t (M GUT ) ~ 1, which reads a 3 (M z ) ~ 0.122, is modified 
to a 3 (M z ) — 0.124 if we keep the same value of T S usy and sin 2 Q w {Mz) but we take 
a smaller value of the top quark Yukawa coupling Y^Mqut) — 0.1. The latter result 
is in agreement with the value quoted in Ref. for the same values of sin 2 9w(Mz) and 
Tsusy- The above analysis also explains the small numerical difference observed in Tables 
5. a and 5.b - sin 2 9 w {Mz) = 0.2335 and 0.2324 - between the lower bounds on Tsusy for 
mb(Mft) = 4.3 and 4.6 GeV, which come in both cases from the upper experimental 
constraint on (Mz) but involving different values of Y^Mqut)- 

In the above discussion, we have assumed that the only splitting in the sparticle 
spectrum is due to the properties of the SUSY particles under SU(3) C . However, as we 
mentioned before, more complicated situations may naturally occur. Nevertheless, since 
for a given value of mft(Mft) and sin 2 9w both Y^Mgut) and a 3 {Mz) are only dependent on 
Tsusy, our results have a more general validity, independent of the exact supersymmetric 
spectrum splitting. 

6 Top Quark and Higgs Mass Predictions 

For the range of values of the running bottom quark mass considered in the present work, 
the requirement of gauge and Yukawa coupling unification implies the convergence of the 
top quark Yukawa coupling to its infrared quasi fixed point value. This yields accurate 
predictions for the running top quark mass, which, however, depend strongly on the range 
of values we take for tan (3. In minimal supergravity models, in which all scalar sparticles 
are assumed to acquire the same soft supersymmetry breaking mass, mo, at high energy 
scales, constraints on tan (3 may be obtained by requiring a proper radiative breaking of 
the SU(2)l x U(l)y symmetry. In these models the characteristic values of the ratio 
of the Higgs vacuum expectation values are in the range 1 < tan/5 < tan/3 M , where 
tan/3 M ~ 30(45) for a top quark mass m t ~ 120(180) GeV, respectively. In the following, 
however, we shall be less restrictive and allow the upper bound on tan f3 to be given by 
the breakdown of the perturbative consistency of the theory in the bottom quark Yukawa 
sector, as discussed in section 4. 
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In Fig. 2 we plot the values of the running top quark mass mt(m t ) as a function of 
tan/3 for different values of m^Mj,) and taking the mean values of sin 2 9w(Mz) = 0.2324 
and a^{Mz) = 0.12. From the above figure it follows that, the smaller the value of the 
bottom quark mass is, the larger the value of m t results for any fixed value of tan f3, in 



agreement with previous one loop results [11]. This behavior is easily understood due to 
the interrelation between the top and bottom quark Yukawa couplings. In addition, the 
existence of the IR quasi fixed point for h t explains the mildness in the actual variation 
of the running top quark mass values. 

In Fig. 3, we present the values of the running top quark mass as a function of tan f3 for 
the considered mean values of sin 2 9w(Mz) and mb(Mb) and for different values of a 3 (M^), 
thus, for different values of Tsusy and Y^Mqut)- For any given value of tan/3, a larger 
value of as(Mz) corresponds to a larger value of mt- Such behavior is predictable due to 
the interrelation between the strong gauge coupling constant and the top quark Yukawa 
coupling. However, it is worth noticing that, although, for example, for tan/3 = 10, a 
variation in the low energy value of the strong coupling constant, a^iMz) — 0.115 — > 0.130 
induces a large variation in Y^Mqut) = 0.28 — > 1, due to the presence of the infrared 
quasi fixed point this translates in an increase in m t of about a five percent of its original 
value. 

Quite generally, for values of tan/3 < 10, the top quark mass is a monotonically 
increasing function of tan /?, varying from m t = for tan j3 = to a value of m t far above 
its lower experimental bound, for tan/3 ~ 10. Hence, the experimental lower bound 
on the top quark mass, allows to put a lower bound on tan/3, below which the model 
is inconsistent with experimental data. As a matter of fact, for the range of values of 
as(Mz) and rrib(Mb) considered in this work, it follows that tan/3 > 0.4. Observe that 
a top quark mass M t ~ 138 GeV, which provides the best fit for the experimental data, 
may only be achieved either for values of tan/? close to one, or for very large values of 
tan p, close to the upper bound on this quantity. 

A phenomenologically interesting possibility is the unification of the three Yukawa 
couplings of the third generation. Since the RG equations of the top and bottom quark 



28 



Yukawa couplings differ only in electroweak terms, the unification of these two couplings 
may only be achieved for values of tan/3 ~ m t (m t ) / 'm b (mt) ■ From the behavior of m t (m t ) 
as a function of tan/? and the fact that m&(m*) ~ 3 GeV, it is easy to see that the 
unification of the three Yukawa couplings occurs for values of tan j3 ~ 60. For example, for 
sin 2 6 W {M Z ) ~ 0.2324, a 3 {M z ) ~ 0.122 and m b (M b ) = 4.1,4.3,4.6 GeV, the unification 
of the tau, bottom and top Yukawa couplings occurs for tan/3 ~ 64,60,53, respectively 
(see tables l.a - l.c). 

Concerning the Higgs sector, as we mentioned in section 5, the light Higgs mass 
predictions are not only dependent on the value of the top quark Yukawa coupling, but 
on the absolute value of the squark mass as well. In fact, for rriq ^> m t , the dominant 
radiative corrections to the tree level Higgs mass value, m| = Mf cos 2/3, are proportional 
to the logarithm of the ratio of the top and stop masses, with a factor which depends on 



the fourth power of the top quark Yukawa coupling |21fl-|23 



m\(m t ) ~ M\ cos 2 2/3 + -^v 2 sin 4 (3 h 4 t (m t ) In ( ^ ). (35) 

2tt 2 \m t J 

Since the unification condition only fixes the value of the effective supersymmetric thresh- 
old scale Tsusyi which is only weakly dependent on the squark masses, even if we fix the 
value of tan /3, the lightest CP even Higgs mass can not be predicted within this scheme. 
However, an upper bound on rriq may always be obtained by the additional requirement 
of the stability of the hierarchy of scales under radiative corrections, or equivalently, the 
absence of fine tuning of the soft supersymmetry breaking parameters in order to achieve 
the proper weak scale. This upper bound is, quite generally, of the order of a few TeV. 
Hence, even though the Higgs mass scale may not be directly determined from the uni- 
fication condition, an upper bound on it may be obtained, which, for a given Tgc/sy is 
achieved for the maximum allowed value of rriq. 

In Figure 4, we plot the Higgs mass as a function of tan/3 for sin 2 9w(M z ) = 0.2324, 
rrib{M b ) ~ 4.3 GeV, T S usy — ITeV, thus a 3 (M z ) ^ 0.115, and for different values of 
the colored and uncolored particles. Observe that, as expected, the Higgs mass increases 
with the colored sparticle mass. For rriq = 1, 3, 10 TeV, the lightest CP even Higgs mass, 
as a function of tan/3 reads = 54, 67, 78 GeV for tan/3 = 1; rrih = 111, 128, 143 GeV 
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for tan/5 = 3 and = 126, 142, 158 GeV for tan (3 = 10, respectively. 

7 Conclusions 

We have analyzed in detail the conditions under which gauge and Yukawa coupling unifi- 
cation take place within the MSSM. We derived an analytical expression for the effective 
super symmetric threshold scale as a function of the supersymmetric particle masses. We 
obtain that the effective scale Tsusy is only slightly dependent on the squark and slepton 
masses, while it is strongly dependent on the Higgsino and gaugino masses. In particular, 
we have shown that, in a first approximation - whenever no large mixing between the 
Higgsinos and the weak gauginos occurs, the gaugino masses proceed from a common soft 
supersymmetry breaking term at Mqut, and all sparticle masses are above Mz - the su- 
persymmetric threshold scale is given by the overall Higgsino mass multiplied by a factor 
which depends on the ratio of the strong and the weak coupling constants. This result is 
quite interesting, due to the fact that, if there is no significant mixing, the Higgsino mass 
is governed by the supersymmetric mass parameter ji appearing in the superpotential and, 
hence, there is no explicit dependence on the soft supersymmetry breaking terms in the 
expression of Tsusy- Quite generally, our analysis is made under the assumption that all 
SUSY particles have masses above M z . We have, however, analyzed the modification of 
the expression for the effective supersymmetric threshold scale in the case in which light 
sparticles appear in the spectrum as well. 

We show that, not only the value of the strong gauge coupling, but also the top quark 
Yukawa sector depends dominantly on Tsusy- For a given set of values for sin 2 9w and the 
bottom quark mass, we derive lower bounds on Tsusy coming either from the experimental 
upper bounds on 0:3 (Mz) or from the requirement of perturbative consistency of the top 
Yukawa sector of the theory. Furthermore, within a simplified scheme, in which we assign 
a common mass m wk to the uncolored sparticles and a mass m col to the colored ones, a 
lower bound on T S usy determines an upper bound on the value of m coh for any given 
value of mis- 
performing a two loop renormalization group analysis, we have also computed the 
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predictions for the top quark and Higgs masses as a function of tan (3 in the framework of 
gauge and bottom-tau Yukawa coupling unification. Quite generally, for different values of 
sin 2 9w(Mz) and rrib(Mb) and for various values of as(Mz), we obtain that the running top 
quark mass has to be below 200 GeV, while a value of m t ~ 140 GeV is only compatible 
with values of tan/3 quite close to one, or with extremely large values of tan/3 ~ 60. 
Contrary to what happens with a^(Mz) and m t , for fixed values of the bottom quark 
mass and the weak mixing angle, the light Higgs mass depends not only on Tsusy, but on 
the value of the top squark mass as well. For ~ 1 TeV, for example, the Higgs mass 
varies from values which are close to its experimental lower bound if tan /3 ~ 1 to values 
of rrih — 130 GeV for tan/3 ~ 10. 

Finally, we should mention that in the present work, we have not considered the effects 
of possible threshold corrections to the gauge and Yukawa couplings arising at the grand 
unification scale. These corrections can only be computed once we know the specific grand 
unified scenario at the high energy scales. In spite of that, under general assumptions, 
a somewhat general parametrization of such corrections may be done. In Ref. 0, for 
example, the high energy scale threshold corrections associated with the gauge coupling 
constants were considered. We should stress that, within the framework of gauge and 
Yukawa coupling unification, not only the threshold corrections to the gauge couplings, 
but also those to the Yukawa couplings must be considered. In fact, if the Yukawa thresh- 
old corrections were as large as 5 to 10%, they would have relevant effects on the top quark 
mass predictions. We shall concentrate on this subject in a forthcoming publication ||19|| . 
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Table Captions. 



Table l.a. Comparison between the predictions obtained within the two loop RG anal- 
ysis, in the framework of gauge and Yukawa coupling unification, and the approximate 
results derived from the infrared quasi fixed point predictions, for sin 2 9w(M z ) = 0.2324, 
a 3 (M z ) ~ 0.122 and m b {M b ) = 4.6 GeV. 

Table l.b. The same as in Table l.a but for m b (M b ) = 4.3 GeV. 
Table I.e. The same as in Table l.a but for m b (M b ) = 4.1 GeV. 

Table 2. Dependence of a 3 (M z ) on sin 2 9 W (M Z ) and on the effective supersymmetric 
threshold scale T S usy in the framework of the two loop RG analysis with gauge and 
Yukawa coupling unification, for m b (M b ) =4.3 GeV. 

Table 3. Dependence of Y^Mgut) on a 3 (M z ) and Tsusy, in the framework of gauge and 
Yukawa coupling unification. The results are for sin 2 9 W (M Z ) = 0.2324 and m b (M b ) = 4.3 
GeV, while fixing m co i = 1 TeV and varying m wk to get the corresponding value of the 
supersymmetric threshold scale. 

Table 4. Top quark Yukawa coupling and Higgs mass dependence on the precise values 
of m co i and m w k for a fixed value of Tsusy = 1 TeV. The results are for sin 2 9w(M z ) = 
0.2324, m b (M b ) = 4.3 GeV and fixing tan (3 = 4. 

Table 5. a. Lower bounds on the effective supersymmetric threshold scale as a function of 
the running bottom quark mass for: sin 2 9w(M z ) = 0.2335. The values of a^{M z ) and 
Yt(McuT) obtained for each particular set of low energy parameters are also given. 

Table 5.b. The same as in Table 5.a. but for sin 2 ^(M z ) = 0.2324. 
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Table 5.c. The same as in Table 5. a. but for sin 2 Ow(Mz) = 0.2315. 
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Figure Captions. 



Fig. l.a. Lower bound on Tsusy as a function of the colored and uncolored sparticle 
mass scales, m co i and m w k, respectively, for a running bottom quark mass m b (M b ) = 4.6 
GeV and for sin 2 9 W (M Z ) = 0.2315 (dashed line), sin 2 6 W (M Z ) = 0.2324 (dot-dashed line) 
and sin 2 6 W (M Z ) = 0.2335 (solid line). 

Fig. l.b. The same as figure l.a. but for mb(M b ) = 4.3 GeV. 
Fig. I.e. The same as figure l.a. but for m b (M b ) = 4.1 GeV. 

Fig. 2. The running top quark mass as a function of tan (3 for sin 2 6 W (M Z ) = 0.2324, 
a 3 (M z ) c± 0.122 and varying the bottom quark mass to be m b (M b ) =4.1 GeV (dot-dashed 
line), m b (M b ) = 4.3 GeV (solid line) and m b (M b ) = 4.6 GeV (dashed line). 

Fig. 3. The running top quark mass as a function of tan/3 for m b (M b ) = 4.3 GeV 
and sin 2 9 W (M Z ) = 0.2324 and varying the value of the strong gauge coupling constant 
to be a 3 (M z ) = 0.115 (solid line), 0.122 (dashed line) and 0.13 (dot-dashed line). 

Fig. 4. The light Higgs mass as a function of tan f3 for sin 2 9w = 0.2324 and m b (M b ) = 4.3 
GeV while varying the value of the colored sparticle mass scale to be m co \ = 1 TeV (dashed 
line), 3 TeV (dot-dashed line) and 10 TeV (solid line) and taking m w k accordingly to fix 
Tsusy = 1 TeV. 
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IR Fixed Point 


RG solution 


rat (tan /3 = 1) [GeV] 


144 


141 


ra t (tan (3 = 3) [GeV] 


192 


184 


ra t (tan /3 = 10) [GeV] 


201 


192 


ra t (/i t = M [GeV] 


188 


170 


tan/9(/it = hb) 


58 


53 



Table l.a 





IR Fixed Point 


RG solution 


ra t (tan /3 = 1) [GeV] 


144 


144 


ra t (tan /3 = 3) [GeV] 


192 


187 


ra t (tan /3 = 10) [GeV] 


201 


196 


ra t (/i t = M [GeV] 


188 


180 


tanfl(h t = hb) 


63 


60 



Table l.b 





IR Fixed Point 


RG solution 


ra t (tan /3 = 1) [GeV] 


144 


144 


ra t (tan (3 = 3) [GeV] 


192 


188 


m t (tan/3 = 10) [GeV] 


201 


197 


m t (h t = h b ) [GeV] 


188 


185 


tan/9(/it = hb) 


66 


64 



Table l.c 



sin 2 9 W (M Z ) 


a 3 (M z ) for Tsusy = 1 TeV 


a 3 (M z ) for T SUSY = 100 GeV 


0.2335 


0.111 


0.118 


0.2324 


0.115 


0.122 


0.2315 


0.118 


0.126 



Table 2. 



m wk [GcV] 


TsusrlGeV] 


a 3 (M z ) 


Y t (M GUT ) 


10 3 


10 3 


0.115 


0.2 


200 


15 


0.127 


0.7 


150 


7 


0.130 


1.0 



Table 3. 



mcoi [TeV] 


m wk [TeV] 


Y t {M GUT ) 


a 3 {M z ) 


m [GeV] 


m h [GeV] 


1 


1 


0.28 


0.115 


188 


118 


3 


1.93 


0.29 


0.115 


191 


135 


10 


3.89 


0.33 


0.115 


193 


150 



Table 4. 



sin 2 6 w (M z ) = 0.2335 


T™ Y [GeV] 


a 3 (M z ) 


Y t (M GUT ) 


m b (M b ) = 4.6 GeV 


2.5 


0.13 


0.3 


m b (M b ) = 4.3 GeV 


1.5 


0.13 


0.9 


m b (M b ) = 4.1 GeV 


12 


0.123 


1.0 



Table 5.a 



sin 2 9 W (M Z ) = 0.2324 


T™ Y [GeV] 


a 3 (M z ) 


Y t (M GUT ) 


m b (M b ) = 4.6 GeV 


9 


0.13 


0.34 


m b (M b ) = 4.3 GeV 


6.5 


0.13 


1.0 


m b (M b ) = 4.1 GeV 


90 


0.122 


1.0 



Table 5.b 
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sm 2 9 w (M z ) = 0.2315 




a 3 (M z ) 


Y t (M GUT ) 


m b (M b ) = 4.6 GcV 


25 


0.13 


0.4 


m b {M b ) = 4.3 GcV 


45 


0.128 


1.0 


m b (M b ) = 4.1 GeV 


350 


0.121 


1.0 



Table 5.c 
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